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Abstract 


The  spectral  analysis  of  a  series  of  equally  spaced  samples  of  a  time¬ 
stationary  process  becomes  difficult  when  samples  are  missing  or  sizable  data 
gaps  occur  within  the  interval  of  interest.  A  linear  prediction  algorithm  can 
be  used  to  fill  in  the  missing  data  with  estimates  that  are  spectrally  con¬ 
sistent  with  the  data  that  are  observed.  Simulated  and  practical  radar 
examples  demonstrate  an  improvement  in  resolution  and  a  reduction  of  sidelobe 
interference  levels.  Computer  programs  are  provided  which  accomplish  the 
extrapolation  and  interpolation  for  complex  data. 
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A. 


Problem  Definition 


When  a  spectral  transformation  of  a  sampled  process  is  performed,  one 
must  account  for  any  samples  that  are  missing.  Assigning  a  value  of  zero  to 
missing  data  prior  to  Fourier  transformation,  for  example,  introduces  false 
frequencies  and  greatly  increases  sidelobe  levels.  Clearly,  an  interpolation 
scheme  is  needed  that  can  cope  with  missing  data  and,  at  the  same  time,  will 
not  degrade  the  spectral  information  contained  in  the  data  that  are  observed. 

Occasional  missing  samples,  well  separated  from  each  other,  can  be 
estimated  with  simple  interpolation  procedures  (polynomial  or  parabolic  fits, 
spline  fits,  etc.).  However,  data  may  be  missing  in  such  quantity  that  con¬ 
ventional  interpolation  is  inadequate;  data  gaps  longer  than  the  periods  of 
the  sinusoidal  components  in  the  data  cannot  be  easily  bridged  with  simple 
functions.  A  more  sophisticated  approach  becomes  necessary,  and  the  use  of  a 
data-adaptive  linear  prediction  filter  is  one  feasible  alternative. 

In  radar  data  processing,  missing  data  or  data  gaps  may  occur  for  a 
variety  of  reasons: 

(a)  hardware  fails  to  transmit  pulses  or  receive  echoes  properly; 

(b)  radar  transmits  when  it  should  be  receiving  echoes  (range 
eclipsing) ; 

(c)  resources  are  saturated  by  many  targets  that  must  be  watched 
simultaneously  (panic) ; 

(d)  burst  waveforms  are  purposely  silent  between  bursts; 

(e)  poor  signal- tc-noise  makes  detections  sporadically  unreliable. 


In  any  case,  the  missing  samples  (in  these  examples,  complex  samples  with 
amplitude  and  phase)  must  be  filled  in  before  Doppler  processing  can  be 
accomplished. 

B.  Description  of  the  Method 

The  use  of  a  linear  prediction  filter  to  extend  a  finite  complex  data 
set  before  Fourier  transformation  was  first  proposed  and  described  by  Bowlinq 
(1977).  Applying  this  original  algorithm,  Tomlinson  and  Ackerson  (1978) 
demonstrated  clutter  and  sidelobe  reduction  in  the  Doppler  processing  of  a 
train  of  radar  pulses. 

In  the  application  of  interest  here,  the  prediction  algorithm  is  used  to 
predict  estimates  of  missing  data  by  extrapolating  from  observed  data.  For 
example,  suppose  an  observation  interval  contains  randomly  missing  samples  and 
gaps.  The  procedure  is  as  follows: 

(1)  Locate  and  designate  the  missing  samples  to  be  estimated. 

(2)  Find  the  longest  continuous  span  of  data  within  which 
there  are  no  missing  samples. 

(3)  Calculate  an  N-point  linear  prediction  filter  from  the 
longest  continuous  span  of  data  found  in  step  (2). 

(4)  Calculate  an  estimate  of  each  missing  sample  immediately 
to  the  left  and  to  the  right  of  the  longest  continuous 
span  of  data  (a  total  of  two  estimates,  one  on  each  side) . 

(5)  Return  to  step  (2)  until  all  missing  data  have  been  esti¬ 
mated.  Note  that  estimates  from  step  (4)  are  to  be  treated 
as  observations  on  an  equal  basis  with  the  original  data. 
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That  is,  the  longest  continuous  span  of  data  is  increasing 
in  length  as  estimates  fill  in  the  holes,  one  by  one,  to 


the  left  and  to  the  right. 

When  the  longest  continous  span  of  data  finally  terminates  at  one  of  the 
endpoints  of  the  observation  interval,  estimates  continue  to  be  made  toward 
the  other  endpoint  until  all  missing  points  have  been  filled  in.  The  length 
of  the  prediction  filter  may  remain  a  constant,  or  vary  according  to  the 
current  length  of  the  longest  span  of  data. 

C.  Simulated  Examples 

A  simple  example  shows  the  improvement  in  the  power  spectrum  of  a  data 
set  containing  missing  samples  and  gaps. 

The  real  and  imaginary  parts  of  a  sampled  sum  of  three  complex  sinusoids 

are  shown  in  Figs.  l(a,b).  No  noise  has  been  added  and  no  samples  are  missing. 
Figure  1(c)  is  the  true  power  spectrum  calculated  with  a  standard  FFT.  No 
weighting  function  has  been  used. 

Now  if  samples  are  randomly  zeroed  out  and  data  gaps  are  introduced  as 
shown  in  Figs.  2(a,b),  the  power  spectrum  in  Fig.  2(c)  shows  increased  side- 
lobe  levels  and  false  frequencies,  both  caused  by  processing  without  estimating 
the  missing  data. 

Figures  3(a,b)  show  the  data  set  after  the  linear  prediction  algorithm  is 
applied,  with  the  power  spectrum  shown  in  Fig.  3(c).  Not  only  do  Figs.  l(a,b) 
overlay  with  3(a,b)  almost  exactly,  but  their  respective  power  spectra  are 
indistinguishable. 
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Another  simple  example  demonstrates  the  performance  of  the  linear  pre¬ 
diction  algorithm  when  data  gaps  occur  periodically,  such  as  is  the  case  for  a 
radar  burst  waveform. 

Figures  4(a,b)  represent  the  process  of  Figs.  l(a,b)  for  which  three  data 
gaps  are  present.  Indeed,  half  of  the  data  are  missing  from  the  observation 
interval,  and  the  gaps  are  longer  than  any  period  exhibited  in  the  data.  The 
power  spectrum  of  Fig.  4(c)  is  a  very  poor  estimate  of  the  true  spectrum 
(Fig.  1(c))  because  no  gaps  have  been  filled  in.  Transforming  only  one  of  the 
short  spans  of  observed  data  gives  a  power  spectrum  with  limited  resolution, 
as  shown  in  Fig.  4(d). 

However,  upon  using  the  prediction  algorithm  on  Figs.  4(a,b),  we  obtain 
Figs.  5(a,b)  and  the  power  spectrum  in  Fig.  5(c),  which  is  an  excellent  esti¬ 
mate  of  the  true  spectrum. 

In  this  case,  the  prediction  algorithm  has  acted  as  a  synergistic  device 
that,  by  linking  short  pieces  of  data  together  with  spectrally  consistent 
estimates,  allows  a  spectral  transform  to  be  performed  over  an  effectively 
longer  piece  of  data.  The  whole,  then,  has  more  resolving  power  than  any  of 
its  parts. 

It  should  be  pointed  out  that  the  data  gaps  need  not  be  periodic  or  equal 
in  length  in  order  for  the  prediction  algorithm  to  fill  them  in. 

D.  Radar  Example 

Radar  is  often  used  to  identify  targets  from  the  time  history  of  the 
velocity  spectrum  of  the  target's  motion  about  its  center  of  mass.  A  series 
or  burst  of  radar  pulses  is  Fourier  analyzed,  and  the  target's  velocity 
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spectrum  is  observed.  If  not  accounted  for  in  the  processing,  missing  pulses 
can  introduce  false  velocity  components  and  lead  to  an  incorrect  character¬ 
ization  of  the  target. 

For  example.  Fig.  6(a)  shows  the  evolution  of  the  velocity  spectrum  of  a 
tumbling  object  for  which  missing  data  and  data  gaps  exist  and  are  set  to  zero 
in  the  radar  pulse  train.  No  estimation  for  the  missing  pulses  has  been  done. 
It  is  therefore  not  clear  if  the  velocities  indicated  are  actually  from  the 
target  or  are  an  artifact  of  the  missing  data.  Figure  6(b)  shows  the  evolu¬ 
tion  of  the  same  velocity  spectrum  upon  using  the  prediction  algorithm  before 
Fourier  transformation.  The  disappearance  of  some  of  the  velocities  cleans  up 
the  spectral  history  and  indicates  which  velocity  components  actually  charac¬ 
terize  the  target. 

E.  Limitations  of  the  Method 

Implicit  in  the  use  of  a  linear  prediction  filter  is  the  assumption  that 
the  data  from  which  the  filter  is  derived  are  coherent.  The  process  being 
sampled  must  be  approximately  stationary  during  the  observation  interval  which 
is  being  analyzed  and  within  which  the  missing  data  and  data  gaps  may  occur. 

Also,  the  prediction  filter  works  best  when  the  spectral  components  are 
approximately  pure  tones,  confined  to  locally  narrow  bandwidths  spaced  within 
the  Nyquist  bounds  of  the  spectral  transform  domain. 

F.  Computer  Codes 

Appendix  I  lists  three  self-contained  subroutines  that  accomplish  the 
method  described  in  Section  B.  GAPFIL,  the  driving  subroutine,  calls  COEFF 
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and  LNPRED,  which  calculate  orediction  filter  coefficients  and  perform  linear 
predictions,  respectively.  Although  the  code  is  desiqned  for  complex  data, 
purely  real  data  can  be  treated  by  setting  the  imaginary  components  to  zero 
before  calling  GAPFIL.  See  Bowling  (1977)  for  a  more  detailed  description 
of  the  linear  prediction  procedure. 

G.  Summary 

This  paper  proposes  the  use  of  a  linear  prediction  algorithm  to  fill  in 
missing  data  and  data  gaps  that  may  occur  within  an  observation  interval  over 
which  a  spectral  transform  is  to  be  made.  False  frequencies  and  sidelobe 
interference,  which  are  artifacts  of  the  missing  samples,  can  be  surpressed  or 
eliminated  by  replacing  the  missing  samples  with  estimates  that  are  spectrally 
consistent  with  neighboring  observed  data.  Large  gaps  can  be  smoothly  bridged 
that  otherwise  could  not  be  satisfactorily  interpolated  by  simpler  schemes. 
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FREQUENCY 


Fig.  1.  (a)  Real  part  of  the  sum  of  three  complex  sinusoids;  no  samples 

are  missing.  (b)  Imaginary  part  of  the  sum  of  three  complex  sinusoids; 
no  samples  are  missing.  (c)  Power  spectrum  of  l(a,b). 
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Fig.  2.  (a)  Real  part  of  Fig.  1  with  randomly  missing  data  and  data  gaps. 

(b)  Imaginary  part  of  Fig.  1  with  randomly  missing  data  and  data  gaps. 

(c)  Power  spectrum  of  2(a,b). 
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Fig.  3.  (a)  Real  part  of  Fig.  2  after  application  of  the  linear  prediction 

algorithm.  (b)  Imaginary  part  of  Fig.  2  after  application  of  the  prediction 
algorithm.  (c)  Power  spectrum  of  3(a,b). 
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FREQUENCY 


Fig.  4.  (a)  Real  part  of  FJg.  1  with  data  gaps.  (b)  Imaginary  part  of  Fig.  1 

with  data  gaps.  (c)  Power  spectrum  of  4(a,b).  (d)  Power  spectrum  of  the 

central  short  data  span  in  4(a,b). 


Fig.  6.  (a)  Doppler  history  of  tumbling  object  when  missing  data  and  data 

gaps  are  not  accounted  for.  (b)  Doppler  history  after  missing  data  are 
filled  in  with  the  linear  prediction  algorithm. 
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APPENDIX  COMPUTER  CODES 

Listings  for  subroutines: 

1)  GAPFIL  -  driving  program 

2)  COEFF  -  filter  coefficient  calculations 

3)  LNPRED  -  estimation  by  linear  prediction 
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n  n 


SUBROUTINE  GAPF1L  (XR  E AL,  XIHAG,  NPTS,  X,  A,  PH,  AA ,  B1,  B2) 


COMPLEX  X,  A,  AA,  B1,  B2,  ARRAY 
DIMENSION  X(1),  A  ( 1)  ,  XHEAL(I),  XIHA3(1) 

DIMENSION  PM  ( 1 )  ,  AA(1),  B1(1).  B2(1) 

DIMENSION  APR  AY  (2  56) ,  HISS(128) 

C 

C 

C  INPUTS  ARE:  IREAL, XIHAG, NPTS 

C 

C  OUTPUTS  ARE  XR E AL , X IH AG , X , A 

C 

c 

C  THE  ARGUMENTS  OP  THIS  SUBROUTINE  APE: 

C 

C  XREAL (NPTS) =  ARRAY  CONTAINING  THE  REAL  PART  OF  THE  DATA 

C  X I HAG ( N PTS) =  ARRAY  CONTAINING  THE  IMAGINARY  PART  OF  THE  DATA 

C  N PTS=  THE  NUMBER  OF  DATA  POINTS  (INCLUDING  HISSING  POINTS) 

C  T  ( NPTS) =  COMPLEX  ARRAY  OF  DATA  POINTS  USED  IN  PREDICTION 

C  A  (NCOEFF) =  COMPLEX  ARRAY  OF  PREDICTION  COEFFICIENTS  GENERATED 

C  WITHIN  THIS  SUBROUTINE  (NCOEFF  BILL  NOT  EXCEED  NPTS/4) 

C  P H ( NCOEFF) *  ARRAY  OF  ERROR  POBER  HISTORY  AS  FILTEP  IS  BUILT 

C  AA  (NCOEFF)  , B1  (NPTS) , B2  (NPTS) =  WORK  ARRAYS 

C 

C  IN  THE  ABOVE  DESCRIPTION,  THE  MINIMUM  SIZES  OF  THE  ARRAYS 

C  AS  USED  IN  THIS  SUBROUTINE  ARE  SPECIFIED 

C 
C 

C  •’'HE  LONGEST  SPAN  OP  CONTIGUOUS  DATA  IS  FIRST  LOCATED  WITHIN  THE 
C  INTERVAL  NPTS.  THEN  THE  HISSING  DATA  IN  THE  ADJACENT  SAMPLE 

C  POSITIONS  AFE  PREDICTED  USING  A  LINEAR  PREDICTION  ALGORITHM. 

C  AS  HISSING  POINTS  AND  GAPS  ARE  FILLED  IN,  THE  BASIS  FROM  WHICH 

C  THE  PREDICTIONS  ARE  BEING  HADE  INCREASES  IN  LENGTH;  NCOEFF 

C  THE  LENGTH  OF  THE  PREDICTION  EILTER,  IS  HOT  ALLOWED  TO  EXCEED 

C  NPTS/4.  THE  DATA  ARRAYS  ARB  RETURNED  WITH  THE  GAPS  AND 

C  HISSING  POINTS  FILLED  IN. 

C 

C 

c 

C  SEE  LINCOLN  REPORT  RHP- 122  FOP  DESCRIPTION  OF  LINEAR  PREDICTION 

C 

C  IF  NPTS  IS  GREATER  THAN  256  ,  INCREASE  DT  HENS  ION  OF 

C  'ARRAY'  TO  ACCOMODATE 

C  IF  MORE  THAN  128  POINTS  ARE  HISSING,  INCREASE  DIMENSION  3F  'HISS'. 

C 

c - 

C 

C  STEP  1  :  IDENTIFY  HISSING  DATA. 

C 

C  HISSING  DATA  APE  THOSE  CONSIDERED  TO  HAVE  VALUES  OF  ZERO.  IF 
C  THIS  CRITERION  IS  NOT  APPROPRIATE  FOR  YOUP  APPIICATION,  THEN 
C  YO'J  WILL  HAVE  TO  REWRITE  ST En  1.  NOHISS  IS  THE  NUMRFP  OF 
C  HISSING  DATA  POINTS,  AND  MISS(J)  IS  THE  LOCATION  OF  THE  JTM 
C  HISSING  POINT.  THUS,  HISS  (J)  CAN  HAVE  ANY  VALUE  BETWEEN  1  AND  NPTS. 
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c 

1*0 
J  =  0 
9  I  *  I*  1 

X ( I)  *  CHPLI (I REAL (I) ,XIHAG(I)) 

IF  (XREAL(I). EQ.O. 0. AID. XIHAG  (I)  .  EQ.O.  0)  J=J*1 
IF  (XRBAL(I). BQ. 0.0. AND. XIHAG(I).  EQ.O. 0)  0ISS(J)*I 
IF  (I. IT. NPTS)  60  TO  9 
NOHISS  *  J 

IP  (NOHISS. EQ.O)  RETURN 

C 

C - - - 

C 

C  STEP  2  :  SEARCH  FOR  LONGEST  SPAN  OF  CONTINUOUS  DATA, 

C  AND  GENERATE  2  NUMBERS  (LI,  L2)  CHARACTERIZING  THE  SPAN. 

C 

LI  =  0 

L2  =  HISS(I) 

HI  =  1 
H2  =  1 

HAX*  HISS  (  1) 

IF  (NOHISS.EQ. 1)  GO  TO  30 

DO  20  J  =2,  NOHISS 

JDIF  =  HISS(J) -HISS  (J- 1) 

IF  (JDIF. LE. H  AX)  GO  TO  20 
HI  =  J-1 
H2  =  J 
L1=HISS(H1) 

L2*NISS  (112) 

HAX*  JDIF 
20  CONTINUE 

30  JDIF  =  NPTS-HISS (NOHISS) 

IF  (JDIF. LE. HAX)  GO  TO  34 
Hi  *  NOHISS 
H2  *  NOHISS* 1 
Li  =  HISS  (HI) 

L2  *  NPTS  ♦  1 
C 

C - 

C 

C  STEP  3  :  DEFINE  NUMBER  OF  LINEAR  PREDICTION  COEFFICIENTS, 

C  AND  PERFORM  PREDICTION. 

C 

34  NCOEFF  *  (L2-L1)/2 

IF  (NCOEFF. LT. 1)  NCOEFF=2 

IF  (NCOEFF. GT.  (NPTS/4)  )  MCGEE''*  (NPTS/4) 

11  *  L  1*  1 

12  *  L2-1 
J  *  0 

DO  40  I  =  II,  12 
J  =  J*  1 

40  AKRAT(J)  *  X(I) 

HP  =  J 
MPP2  =  .1*2 
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n  n  n  o  o  r> 


CALL  COBFP (BP,  ARRAY,  RCOEFP,  A,  PH,  PO,  A  A,  B1,  B 2) 
CALL  LIPRED (H  P,  BPP2,  ARRAY,  HCOIFF,  A) 


STEP  A  :  RB-BSTABLXSH  RBAL  ARD  IHAGIMARY  ARRAYS. 

IF  (L1.GT.O)  XREAL(LI)  -  RBAL  (AIR AY (1 )  ) 

IF  (L1.GT.0)  XIHAG(LI)  -AIBAG (ARRAY (1 ) ) 

IF  (L1.GT.O)  X (LI)  =  ARRAY  (1) 

C 

IF  (L2.LB.VPTS)  XRBAL (L2)  ■  RBAL( ARRAI (HPP2) ) 

IF  (L2.LB.HPTS)  IIHAG(L2)  « AlflAG ( ABRAI (HPP 2) ) 

IF  (L2.LB.HPTS)  X  (L2)  «  ARRAY  (HPP2) 

C 

HI  «  HI-1 

IF  (H1.GE.1)  L  1*HISS  (HI) 

IF  (H1.LT.1)  L1«0 
H2  *  H  2* 1 

IF  (H2  .  LB .  ROHISS)  L2=SISS(H2) 

IF  (H2.GT. ROHISS)  L2«RPTS»1 
IF  (fll.LT.  1.ARD.S2.GT. ROHISS)  GO  TO  SO 
GO  TO  34 

C 

50  RRTORR 
END 
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SUBROUTINE  CO EFP ( NPTS , I , NCOEF P, A, PH, PO , AA , B 1 , B2 ) 


C 

C  THIS  SUBROUTINE  CALCULATES  THE  COMPLEX  BURG  COEFFICIENTS 
C  (HERE  CALLED  ARRAY  A).  THE  ALGORITHMS  USED  HERE  ABE  AR  EXTENSION  OP 
C  ALGORITHMS  DESCRIBED  BY  ANDERSON  (GEOPHY SICS, YOL  39,PEB.  197*) 

C  TO  THE  CASE  OP  A  COMPLEX  SERIES. 

C 

C  PROGRAMMED  BY  S. B.  BONLING,  M IT- LI RCOLH  LABORATORY,  SEPT.  1976. 

C 

COMPLEX  X,A,AA,B1,B2,XBOM,DEM,TMO 
DIMENSION  X(  1)  ,  A(  1)  ,AA  (  1)  ,B1  (  1)  ,B2  (1)  ,  PM  (1) 

TNO=CH  PLI (2. 0,0.0) 

PO=0.0 

DO  10  IT= 1 , NPTS 

DUHHY=  X ( IT) *CO  N JG  (X (IT) ) 

10  PO=  PO*  DUMMY 

PO=  PO/PLOAT (NPTS) 

NM 1 =NPTS- 1 
B1  (1)=X(1) 

B2  (NM1) -X  (NPTS) 

DO  20  IT=2 , NH 1 

B 1  (IT)  =  X  ( IT) 

ITH 1=1 T- 1 

20  B2  (ITH  1)  =X  (IT) 

DO  50  H= 1 , NCOEPP 

HH1»H-1 

NHH=NPTS-M 

IF  (M  .  EQ.  1)  GO  TO  25 
DO  21  IT=  1  ,MM  1 

21  A  A  (IT)  =A  (IT) 

DO  22  IT= 1 , NHM 

B  1  (IT)  =  B  1  (IT)  -COHJG  (A  A  (MM  1)  )  *B2(IT) 

22  B2  (IT) =  B2 (IT*1) -AA  (HHl) *B1  (IT*1) 

25  XNOH=CMPLX (0.0,0. 0) 

DEN=CH  PLX (0. 0,0.0) 

DO  30  I T= 1 , NHM 

X  NOH  =  X  NOM  ♦  B2(IT) *CONJG (B 1 (IT) ) 

30  DEN=DE  N  ♦  B1 ( IT) *CON JG ( B 1 ( IT)  )  ♦  B2  ( I T ) *CO N JG  ( B2  ( IT) ) 

IF  (  REAL  (DEN)  .EQ.  0.0)  GO  TO  35 
A ( N) =  TNO* (XNOH/DEM) 

GO  TO  36 

35  A(H)  *CHPLX(0.  0,0.0) 

36  POREB=PO 

IP  (M  .GT.  1)  PONER=PH  (H-1) 

D0MHY=  A  (M)  *CON JG  ( A  (H) ) 

PH  (M)=PORBR*(1.0  -  DUMMY) 

IP  (H  . EQ.  1)  GO  TO  50 
DO  40  I T*  1 , MM  1 

40  A  (IT) =  AA (IT) -A (H) *CONJG  (AA (M-IT) ) 

50  CONTINUE 
RETURN 
END 
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SUBBOUTIN  E  LB PRED  ( H 1  , N2  , X , NCOEFF, A) 

C 

C  THIS  SUBBOOTIBE  LINEABLT  EXTENDS  T8I  COB PL EX  DkTk  ARRAT  'X* 

C  FROM  HI  POIHTS  TO  H2  POIBTS.  THE  ORIGINAL  DATA  IS  CENTERED 
C  IN  AN  ABRAT  OF  R2  ELEMENTS  AND  BOTH  A  FOBHARD  AND  BACKRARD 
C  EXTENSION  ABE  PEBFOR  MED  UNTIL  THE  TOTAL  BOBBER  OF  POINTS  IS  N2. 
C 

C  N 1=  ORIGINAL  RUBBER  OF  POINTS  IN  ARRAT  X 
C  N2=BUHBER  OF  POINTS  TO  BHICH  ARRAT  X  IS  EXTENDED 
C  X-  COMPLEX  ARRAT  OF  DATA  SAMPLES,  DIMENSION  AT  LEAST  N2 
C  NCOEFF=  NUMBER  OF  PREDICTION  FILTER  COEFFICIENTS 
C  A=  COMPLEX  ARRAT  OF  FILTER  COEFFICIENTS,  DIMENSION  NCOEFF 
C 

C  PROGRAMMED  BT  S. B.  BOWLING,  HIT-LINCOLN  LABORATORT,  SEPT.  1976. 
C 

COMPLEX  X , A 
DIMENSION  X ( 1 )  , A  (1) 

SET  UP  LIMITS  FOB  DO  LOOPS 

L 1=  R2/2  -  N 1/2 

L2  =  N2/2  ♦  N 1/2 

IF  (  HO D (R  1 , 2)  .EQ.  1  )  L2*L2*1 

SHIFT  ORIGINAL  DATA  TO  MIDDLE  OF  ARRAT  X 

DO  100  1=1, HI 
J=  N 1  -  (1-1) 

K=  L2  -  (1-1) 

100  X  (K) =X  (J) 

DO  FORWARD  PREDICTION 

N3=  N2-L2 
DO  200  1=1  ,N3 
J=  L2+ I 

X (J) =C  MPLX (0. 0,0.0) 

DO  200  K=1, NCOEFF 
200  X(J)=  X(J)  ♦  A(K)*X(J-K) 

DO  BACKWARD  PREDICTION 

DO  300  1=1, LI 
J=  LI-  (1-1) 

X  (J) =CHPLX  (0.0,0. 0) 

DO  300  K=1, NCOEFF 

300  X(J)=X(J)  ♦  CONJG  (A  (K)  )  *X  (J*K) 

C 

RETURN 

END 


twv:  s  - 
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